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SUMMARY 

The series describing the perturbations in the orbital 
elements of an artificial satellite due to solar and lunar 
tides are presented in a form convenient for computation. 
In addition, the potential function and its gradient are given 
in terms of rectangular coordinates. 
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LIST OF SYMBOLS 


A - coefficient of trigonometric term in the determining function 

Ap - coefficient of trigonometric term in the disturbing function 

A^q - factor in the coefficient of the trigonometric term of the determining function depending 
only on the inclination of the satellite orbit plane to the equator 

a - semi- major axis of the satellite 

V " coefficient of trigonometric term depending on the semi-major axis, the eccentricity 
of the satellite, and the mean motion of the argument of the trigonometric term 

c - cosine of half the angle of inclination of satellite orbit plane to the equator 

e - eccentricity of the satellite 

G - Delaunay variable conjugate to argument of perigee of satellite 
g - argument of perigee of satellite 

H - Delaunay variable conjugate to right ascension of node of satellite 
h - Delaunay symbol for right ascension of the node - same as a 
I - inclination of satellite orbit plane to the equator 
I ' - inclination of orbit of disturbing body to the equator 
i ' - index associated with mean anomaly of disturbing body 

- second zonal harmonic of earth 

- Love number 

L - Delaunay variable conjugate to mean anomaly of satellite 
^ - mean anomaly of satellite 

- mean anomaly of disturbing body, including phase angle due to tidal friction 

- mean anomaly of the moon 

- mean anomaly of the sun 

m' - ratio of mass of disturbing body to the sum of the mass of the disturbing body and the 
mass of the earth 

n - mean motion of the satellite 

n ' - mean motion of the disturbing body 

q' - index associated with argument of perigee of the disturbing body 
R - mean radius of earth 
Rp - periodic part of the disturbing function 
r - geocentric distance of the satellite 
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**(D ’ **© “ geocentric distance of fictitious moon (sun) 

S - determining function for solar and lunar tides 

^q' i' V " coefficient of trigonometric term in the determining function depending on elements of 
disturbing body only 

a - argument in the disturbing function 

a - mean motion of a 

- geocentric angles between r and , t* 

y - sine of one half the angle of inclination of satellite orbit to the equator 
€ - phase lag due to tides associated with the right ascension 

- phase angle associated with mean anomaly of the moon due to tidal friction 
K 2 ^ phase angle associated with mean anomaly of the sun due to tidal friction 

V - index associated with right ascension of the node 
0 - right ascension of the satellite 

- right ascension of the moon 
n - mean motion of O 

- mean motion of 

- argument of perigee of disturbing body 

0 - difference of right ascension of the satellite and the disturbing body, including the phase 
lag 6 

0 - mean motion of 0 

A6 - product of gravitational constant and mass of the earth 
k 2 - gravitational constant 
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EFFECTS OF THE SOLAR AND 
MOTION OF AN ARTIFICIAL 


LUNAR TIDES ON THE 
EARTH SATELLITE 


INTRODUCTION 

In References 1 and 2, it has been shown that solar and lunar tides may cause perturbations 
in the trajectories of satellites; also, the disturbing potential which describes the variation of the 
earth’s geoid with time has been discussed. 

In this report, a determining function corresponding to the disturbing function for the tides is 
introduced, from which perturbations in the orbital elements are derived. The perturbations are 
given in a form suitable for programming. 

It is also useful to have available the formulation of the potential function and its gradient in 
terms of rectangular coordinates. These are presented in the last section. 


THE DETERMINING FUNCTION 

The periodic part of the disturbing function is given in Reference 1 in the form 


R 


p 



cos a 


( 1 ) 


where Ap is a function of the coordinate elements of the satellite and the disturbing body, and a is 
defined by 


a ~ 'V^ + q'o)' + v6 . (2) 

Here, q', i', and v are integers defining the indices of the summation (1). The index q' takes 
on only two values, zero and two, and thus generates Tables 1 and 2 at the end of the report. When 
q' equals zero, i' takes the values -1, 0, and 1, while v takes on the values 0, 1, and 2. When q' 
equals two, i' assumes the values -1, -2, -3, and -4, while v varies from -2 to +2. The case when 
all three integers are zero simultaneously leads to secular terms and is discussed in Reference 1. 

The angular variables in Equation (2) are defined by 

~ ~ moon 6 ~ Q - ~ € for the moon 

( 3 ) 

* ~ ~ ^2 ^ ~ ^ ^ the sun , 
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where 

= mean anomaly of the moon 
= mean anomaly of the sun 
= right ascension of the node of the moon 
= right ascension of the node of the satellite 
o)' = argument of perigee of the disturbing body 
^ 1 , ^ phase angles due to tidal friction. 

From Equation (1), we find by the methods of Reference 3, for example, that the first order 
determining function S is given by 


S 



A s in a 


(4) 


where 

A = Ap/i 

a = mean motion of a. 

The quantity A in Equation (4) may be decomposed into the factors 


(5) 


where 

kj = Love number (approximately .3) 

R = mean radius of the earth 
n ' = mean motion of the disturbing body 

m' = ratio of the mass of the disturbing body to the sum of the mass of the disturbing body 
and the mass of the earth. 

The function is independent of the elements of the satellite and is given by the formula 


T / 

q 


a;., 


b:- 


( 6 ) 


A' - and B!- - are functions of the elements of the disturbing body only and are given in Table 1 for 

uq 1 q 

q' = 0 and Table 2 for q' = 2. 

Apj is a fimction of y only, where 

y = sin Y (7) 
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and I is the inclination of the orbit plane to the earth's equator. is given in Tables 1 and 2 for 
the appropriate ranges of the index v. can be written explicitly as 




,3(i_e2)3/2^ 


( 8 ) 


where a and e are the semi-major axis and eccentricity of the satellite. The motion a is given by 

a = i'n'+q'^'+v^ , (9) 

where n', w* , and 0 are the mean motions of 'C'*, oj', and ^ defined previously. Also, 


e 

e 


n - for the moon 


O for the sun 


n - 


3 n J2 cos I 
2 a2(l-e2)3 


( 10 ) 


( 11 ) 


where 

J 2 = second zonal harmonic of earth* s potential 
n = mean motion of the satellite. 


THE PERTURBATIONS IN THE ELEMENTS 
If we introduce the Delaunay variables: 

L = (/xa)>/^ 

G = l(i 
H = G cos I 

= mean anomaly of the satellite 
g = argument of perigee of the satellite 
h = n = right ascension of the node of the satellite, 

where m = product of the gravitational constant and the mass of the satellite, we can derive the 
perturbations in the elements from the determining function defined by Equation 4. 
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From the theory of canonical variables, and since the determining function S is independent 
of I and g, the Delaunay variables L and G (and consequently a and e) are not affected by the solar 
or lunar tides. The perturbations in the remaining elements are given by 

8H - n (l - sin I Si = +k2R^n'^m' ^ ^ i'v V ^ 


Sg 


as 

-k,R5n'^m' 


dC • 

q 1 V 
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dL 

dL 
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q 1 

sin a 
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-kjRS n'^m' 

p 

ac 

q i V 

> 
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aG 


c • . • 

q 1 v 

sin a 

as 

-kjRS n'^m' 2^ 

p 

ac /./ 

q 1 V 
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an 


dH 


C / . / 

q 1 V 

sin a 


( 12 ) 


In order to facilitate deriving differential coefficients, it is useful to note that Equation (11) in 
terms of Delaunay variables becomes 


^ 3 

2 


(13) 


From Equations (2), (5), and (13) we find 
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an 
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(14) 


From the definition of given in (8) and the Equations (14), we find 
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q 1 V 
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n a^ (l - cos I a 
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Furthermore, since is a function of y only, we have 


dL 


dG 


dk. 


dH 




dk. 


4n (l - e^) 7 


(16) 


vQ 


dA„ 


4n - e ^ j 7 


where c - cos 1/2 . 

Combining the above results, the perturbations in the elements can be written in the form 


" " Zy Ko f cos a 


S'l = F(l-e^) 


"‘L 




5g ^ F 






2 2 
C-i - 7 ^ 


dk. 


4y ^q'iV dy 


(17) 


8h = F 


I 


V n 


dk. 


-7^ a 47 dy 


where F is the constant introduced in Reference 1 and is given by 


k2 n ' ^ 

n (1 - e^) ^ 


(18) 


THE DISTURBING FUNCTION AND ITS GRADIENT IN RECTANGULAR COORDINATES 
Neglecting the parallactic part, the disturbing function due to tidal effects is 


^ _ k2R5 kj f 


2^3 1 r *3 (3 cos 2 /3 * - 1 ) + *°3 (3 cos ^ /3* - l) 

L c 
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With respect to an orthogonal, earth- centered inertial coordinate system where the z-axis coin- 
cides with the earth’s axis of rotation, let 


(X, y, z) 

(Xc , Vc . Zci 

; *) 
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* _ / * * 
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Table 1 
(q' = 0) 


A'.,' 

0 

dA.o 

dy 

V 

4 (l - 6y' ^ + 67' 

1 “ 67^ + 67"* 

-12y(l -2r^) 

0 

3c' 7' (1 - 27 '^) 

C7 (1 - 27 ^) 

c- ‘ (1 - 8r^ + 87'') 

1 

3c ' 2 7' ^ 

c2 72 

2y(l -2y^) 

2 


BS / 
1 q 


i' 


3 , 
2 ® 


-1 


1 


0 


3 , 
2 ® 


1 


I 


Table 2 


(q' = 2) 
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